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Abstract: Trees growing in a semi-arid sandy environment are often exposed to drought conditions due to 
seasonal variations in precipitation, low soil water retention and deep groundwater level. However, 
adaptability and plasticity of individuals to the changing drought conditions greatly vary among tree 
species. In this study, we estimated water use (T;) of Mongolian Scots pine (MSP; Pinus sylvestris var. 
mongolica Litv.) based on sap flux density measurements over four successive years (2013—2016) that 
exhibited significant fluctuations in precipitation in a semi-arid sandy environment of northern China. 
The results showed that fluctuations in daily T; synchronously varied with dry-wet cycles of soil moisture 
over the study period. The daily ratio of water use to reference evapotranspiration (T;/ETo) on sunny days 
in each year showed a negative linear relationship with the severity of drought in the upper soil layer (0—1 
m; P<0.01). The decrease in Ts induced by erratic drought during the growing season recovered due to 
precipitation. However, this recovery ability failed under prolonged and severe droughts. The T;/ET» ratio 
significantly declined with the progressive reduction in the groundwater level (gw) over the study period 
(P<0.01). We concluded that the upper soil layer contributed the most to the T, of MSP during the 
growing season. The severity and duration of droughts in this layer greatly reduced T;. Nevertheless, gw 
determined whether the T, could completely recover after the alleviation of long-term soil drought. These 
results provide practical information for optimizing MSP management to stop ongoing degradation in the 
semi-arid sandy environments. 
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1 Introduction 


Reforestation has been widely used in semi-arid areas to control soil erosion and to capture 
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carbon for regional sustainability (D'Odorico and Porporato, 2006). However, tree growth in 
water-limited ecosystems is often exposed to significant environment-related challenges, such as 
drought (Mereu et al., 2009; Wesche et al., 2011; Su et al., 2014). Many environmental factors 
influence the survival and growth of trees in semi-arid areas, for instance, amount and variability 
of precipitation, soil water storage and availability of groundwater (Meinzer et al., 2006; Zhang 
Jinhu et al., 2018). Increasing frequency and severity of drought are known to have severe 
negative impacts on forest ecosystems and increase mortality rate of trees (Leo et al., 2013; Barbeta 
et al., 2015; Gazol et al., 2018). Impacts of drought on ecological and hydrological processes 
associated with trees in fragile ecosystems have attracted sustained and intensive attention 
(Bovard et al., 2005). Quantification of tree water use and its long-term response to drought can 
aid in understanding of how trees adapt to a changing environment. 

Mongolian Scots pine (MSP; Pinus sylvestris var. mongolica Litv.) is widely distributed in the 
Da Higgan Ling Mountains (50°10'-53°33'N, 121°11'-127°10'E) of northern China and in parts 
of Russia and Mongolia (Zhu et al., 2008; Zheng et al., 2012). MSP is a popular tree species for 
reforestation in northern China due to its high tolerance to both drought and low temperature. 
Total area of MSP plantation established by the Three-North Shelter Forest Program since 1978 
exceeds 6.7x10° hm? in northern China (Zheng et al., 2012), which greatly contributes to the 
control of desertification. However, severe degradation in the artificial forests of MSP has 
occurred since the mid-1990s, and this degradation is characterized by poor health and a large 
number of deaths of MSP, particularly in the Horqin Sandy Land in northern China (Jiao, 2001; 
Zhu et al., 2008). 

Degradation in the artificial forests of MSP in the Horqin Sandy Land is mainly attributed to 
the high inter-annual and intra-seasonal variations in precipitation and declining groundwater 
level (Jiao, 1989; Song et al., 2014). These impacts seems reasonable because MSP is a 
shallow-rooted species with more than 85% of roots distributed in the upper 0—40 cm soil layer in 
this region, and root density sharply decreased below 100 cm soil layer (Su et al., 2006), which 
limits water absorption from deep soil. However, without frequent and long-term monitoring of 
tree water use and environmental variables, we still have limited knowledge of how and to what 
extent precipitation-dominated water availability is responsible for the degradation of MSP. 
Therefore, this study aims to: (1) quantify sap flux density and daily water use of MSP over four 
successive years (2013-2016) with contrasting precipitation fluctuations; (2) determine the 
linkage between variation in the daily tree water use and intra-seasonal dry-wet cycles of soil 
moisture; and (3) investigate the ability of tree water use recovery after long-term successive and 
severe droughts. 


2 Materials and methods 


2.1 Study area 


The study was carried out at the Zhanggutai National Desertification Control Trial Station in the 
Horgin Sandy Land, northern China (42°43’'N, 122°22'E; 226.5 m a.s.l.). The experiment was 
conducted in the artificial forests (an area of 40 hm?) with 35-year-old MSPs during 2013-2016. 
Tree density was 625 plants/hm? with a 4-m row spacing between individual trees. The study area 
has a semi-arid, continental climate with an annual mean temperature of 7.9°C, a frost-free period 
of 150-160 d, a mean annual pan evaporation of 1553 mm and an annual precipitation of 475 mm 
(with 60%-—70% falling from June to August) over the last 30 a (1983-2012; Zhu et al., 2005). 

After 1983, there were 9 consecutive dry years from 1996 to 2004 with below mean value of 
the annual precipitation (Fig. la). The annual mean temperature over the last 34 years has 
significantly increased at a rate of 0.03°C/a (Fig. 1b). 

Soil has a sedimentary aeolian sand layer that exceeds 3 m and an ancient alluvial sand layer 
that exceeds 126 m (Jiao, 1989). Bulk densities in the 0.0-0.4 and 0.0-2.0 m soil layers are 1.47 
and 1.57 g/cm’, respectively. Soils in the 0.0-0.4 m layer comprise 87% sand (>0.050 mm), 9% 
silt (0.050—0.002 mm) and 4% clay (<0.002 mm). Sand, silt and clay account for 91%, 4% and 


5% in the 0.0-2.0 m soil layer, respectively. Organic matter content in the 0.0-0.4 cm soil layer is 
approximately 1.54 g/kg and decreases to 0.76 g/kg in the 0.0-2.0 m soil layer. The understory 
plant species are Acer pictum subsp. mono Maxim, Crataegus pinnatifida var. major N. E. Brown., 
Lespedeza bicolor Turcz., Artemisia halodendron Turcz et Bess. and Cleistogenes chinensis 
Maxim. 
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Fig. 1 Annual variations in (a) precipitation and (b) air temperature at the Zhanggutai National Desertification 
Control Trial Station in the Horqin Sandy Land, northern China during 1983—2016. Gray columns show the data 


before the experiment (1983-2012) and black columns are the data of the experiment from 2013 to 2016. Dashed 
lines are linear regressions over the entire period. 


2.2 Experimental design and measurements 


A 400 m? plot (20 mx20 m) containing 25 even-aged trees surrounded by a wire fence was 
selected (Fig. 2). Growth of trees in the experiment was normal in 2013 and 2014, but in 2015 and 
2016, leaves from the upper canopy turned slightly gray. Information of sampled trees is given in 
Table 1. Due to the damage of instrument installed in the trunk, the number of tree samples 
decreased to 11 in 2016. Diameter at breast height (DBH) was measured by a diameter tape, and 
tree height was measured by an altimeter after the growing period. Thicknesses of bark, sapwood 
and heartwood were measured by a sampling core near the installed position of sensors with a 
Pressler incremental borer after the experiment. We measured the thickness by a Vernier caliper 
with tissue boundaries identified based on the color that is a yellow-white sapwood tissue and a 
tan heartwood tissue. We used methyl orange solution to help define the interface where boundary 
was indistinct. DBH (P=0.59) and sapwood width (P=0.55) of sampled trees during 2013-2016 
were not significantly different, indicating a satisfactory uniformity of sample trees. 
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Fig. 2 Sketch of sample trees (gray stars) planted in a 4 mx4 m spaced square grid. Locations of soil moisture 
sensors and data logger are shown. Dashed line is the border fence. 
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Table 1 Information of sample trees during 2013-2016 
DBH (cm) 


Year H (m) H (m) SW (cm) A, (cm?) 
Mean+SD Qı Q3 
2013 (n=25) 18.0+2.7° 16.1 18.9 10.3+0.7* 4.5+1.1° 5.5+0.6° 203+58.6° 
2014 (n=22) 17.142.1° 15.6 18.7 9.3+0.8* 3.7+0.5* 5.3+0.5* 182+42.0* 
2015 (n=13) 17.742.2° 17.2 18.8 9.1+0.8" 3.7+0.3°* 5.4+0.5* 189+452.3°* 
2016 (n=11) 17.741.6° 16.8 18.8 9.5+0.5* 4.0+0.1* 5.4+0.5* 194+431.8* 


Note: DBH, diameter at breast height; Q,, 1‘ quartile of DBH; Q3, 3" quartile of DBH; H, tree height; H», height of the first live branch; 
SW, sapwood width; As, sapwood area; SD, standard deviation. n is the number of sample trees. Different lowercase letters within a 
column among different years indicate significant difference at P<0.05 level. 


2.2.1 Meteorological variable measurements 

Meteorological variables, including global solar radiation (Rs), net radiation (Rn), air temperature 
(Ta), relative humidity (RH), wind speed (u) and precipitation (P), were measured using an 
automatic weather station (ARS, Avalon Scientific, Inc., USA) located 50 m away from the 
experimental plot. All sensors were installed at 2.0 m above the ground except the rain gauge that 
was placed at 0.5 m above the ground. Variables were measured at 10-min intervals and then 
averaged and recorded per hour. We calculated reference evapotranspiration (ETo) using the FAO 
Penman—Monteith equation based on the variables of Rn, Ta, RH and u on an hourly data (Eq. 1; 
Allen et al., 1998): 
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where ETo is the reference evapotranspiration (mm/h); A is the slope of saturated water vapor 
pressure against air temperature T, (kPa/°C); Rn is the net radiation (MJ/m7); G is the soil heat 
flux (MJ/m?); y is the psychrometric constant (kPa/°C); es is the saturated vapor pressure (kPa); ea 
is the actual vapor pressure (kPa); and u2 is the mean wind speed at a 2.0-m height (m/s). 

2.2.2 Soil moisture and groundwater level measurements 

Three sets of soil moisture observation device were equipped at three locations in the 
experimental plot (Fig. 2). Soil moisture sensors (ECH20 EC-5 sensor, Decagon Devices Inc., 
USA) of a set were inserted at soil depths of 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.6 and 2.0 m in each 
location. Soil moisture (0, %) was measured at 10-min intervals with hourly means recorded by a 
SQ2020 data logger (Grant Instruments Ltd., UK). We calibrated the data using a site-specific 
equation based on oven-drying method. Relative extractable soil water (REW, dimensionless) was 
calculated using Equation 2 (Granier, 1987): 


REW = L an, (2) 
Ox i Orin 


where Of (%) is the field capacity that was 18.0% in the 0.0-2.0 m soil depth based on field 
observations; Amin is the minimum soil moisture during the experimental period (2.3%); and @ is 
the mean value of soil moisture from the corresponding soil layer (%). We classified soil droughts 
based on REW in the top 1 m soil layer and growth traits of MSPs obtained from other studies 
near the experimental plot, and classifications of soil drought are given in Table 2. We manually 
monitored groundwater level (gw) in situ once a month. 


2.2.3 Sap flow measurement and tree water use calculation 


Sap flux density in the outermost sapwood layer (0-3 cm width, Js-outer, cm/s) was continuously 
measured using Granier-type thermal dissipation probes (Dynamax Inc., Houston, TX, USA). Sap 
flow sensors were installed in each tree in 2013. Each probe was installed under the cambium on 
northern side of the stem at breast height (1.30 m above ground) with pairs of probes spaced 0.04 
m apart vertically. The upper probe was heated and the lower probe was not heated for a reference. 
We carefully removed each sensor at the end of each growing season (November) and reinstalled 


Table 2 Classification of soil drought based on relative extractable soil water (REW) and referenced growth 
trait of Mongolian Scots pine (MSP) 


Degree of drought REW Soil moisture (%) Growth trait of MSP Reference 
Non-drought REW _; m>0.31 Ô o-1m>0.40% Normal growth Jiao (2001) 
Slight drought 0.20<REW 1 m <0.31 0.36:< @ o1 m <0.40 fe Weak growth Zhu et al. (2005) 
Moderate drought 0.08<REWo_; m<0.20 0.20%< O 0-1 m<0.3 Of 30% of leaves withered Zhu et al. (2005) 
Severe drought REW _; m<0.08 a 0-1 m<O.20 5 Withered leaves and some dead Tang et al. (2015) 


branches 


Note: @ 0-1 mis the mean soil moisture in the 0-1 m soil layer (%). Or is the field capacity (%). 


near the initial position at the beginning of the next growing season (April). Temperature 
difference between the upper (heated) probe and the lower (referenced) probe was measured at 
1-min intervals, with mean values recorded at 10-min intervals using SQ2020 data loggers. The 
sensors were shielded with a thick aluminum-faced foam to minimize warming by radiation and 
exposure to rain as well as physical damage. The Granier empirical equation for Js-outer is given in 
Equation 3: 

on 19x19 4( Ato AP 281, (3) 
AT 

where AT is the actual temperature difference observed between heated and referenced probes 
(°C); and ATo is the maximum AT value (°C) when sap flow is close to zero, which was 
determined over approximately 10 consecutive days by linear regression (Lu et al., 2004; Dang et 
al., 2014). 

Sap flux density in the inner part of sapwood (beyond 3 cm width, Js-inner, cm/s) is low due to a 
relative inactivity of xylem, so a coefficient of 0.56 from Scots pine (Pinus sylvestris L.) was 
adopted following Nadezhdina et al. (2002) and Lu et al. (2004) to estimate Js-inner. 

Sap flow (J, cm*/h) of whole tree was calculated by multiplying sap flux density by 
corresponding sapwood area on an hourly scale. Then, hourly sap flow was summed to the daily 
data (Eq. 4). Daily J; was then used to estimate daily tree water use (daily Ts, mm/d) using 
Equation 5. Considering that trees were at the same age and evenly spaced, we estimated daily Ts 
as the product of average sap flux density of all sample trees with cumulative sapwood area 
(As-stand, M°) divided by total stand ground area (Ag-stana, M?) according to Clausnitzer et al. (2011). 
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where Ji; is the hourly sap flow of tree i at hour j in a day (cm*/h); Js-outer,j and Js-inner,ij are the sap 
flux densities in the outermost sapwood layer and inner sapwood layer, respectively, of tree i at 
hour j in a day (cm/h); As-outer,i and As-inner,i represent the sapwood area of outer 3 cm width and the 
remaining sapwood area of tree i, respectively (m7); As; equals the sum of As-outer,i and As-inner,i 
(m°); n is the numbers of sample trees in each year; As-stana is the total sapwood area of the 25 
sample trees (m7); and Ag-stana is the total ground area of the plot (mô). 


2.3 Statistical analyses 


We normalized hourly and daily Ts values by their maximum values over the entire experimental 
period (i.e., Ts/Ts-max) according to Whitley et al. (2008) and Clausnitzer et al. (2011) for further 
comparison of relative tree water use in different years. We calculated the ratio of daily tree water 
use to reference evapotranspiration (T;/ETo) according to Łabędzki and Kanecka (2009) and Fu et 
al. (2016) to eliminate the influence of meteorological factor. Effect of soil moisture on Ts/Ts-max OT 
T;/ETo was tested by one-way analysis of variance (ANOVA) and a Tukey honestly significant 
difference (HSD) test using SPSS 20 (SPSS Inc., Chicago, IL, USA). Significant correlations between 
Ts/Ts-max Or T;/ETo and environmental variables were determined by Pearson's correction 
coefficient tests at P<0.05 or P<0.01 level. OriginPro 2020 version 9.7 (OriginLan Inc., 
Northampton, MA, USA) was used for plotting. 
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3 Results 


3.1 Seasonal dynamics of daily tree water use and environmental variables 


Amounts of precipitation in 2013, 2014, 2015 and 2016 were 554.0, 384.0, 408.0 and 623.6 mm, 
respectively, showing great inter-annual variations. Using mean annual precipitation (475.0 mm) 
as a baseline, we calculated precipitation anomalies over the study period. We defined 2013 as a 
normal year, 2014 and 2015 as dry years, and 2016 as a wet year. However, precipitation in a year 
is concentrated over quite short periods, resulting in high inter-monthly variations in water supply 
(Fig. 3a). Thus, high amount of annual precipitation does not always guarantee enough available 
soil water for plants during the growing season. 

REW exhibited substantial temporal variations. Heterogeneity of REW between the upper 0-1 
m soil layer (REWo-1 m) and lower 1-2 m soil layer (REW1-2 m) was significant (P<0.01), with 
REWo-1 m being higher than REW1-2 m. REW decreased from 2013 to 2015, especially in 2014. 
Values of REWo-1 m in the dry years (2014 and 2015) ranged from 0.29 to 0.30 that is distinctively 
lower than those of in the normal year (2013; 0.48) and the wet year (2016; 0.50). Statistically, 
number of days under drought (i.e., REWo-1 m<0.31; Table 2) in 2013, 2014, 2015 and 2016 
accounted for 8%, 62%, 50% and 16% of the whole year, respectively. In particular, the number 
of severe drought was 19 d, followed by 45 d moderate drought in 2015 (Fig. 3b). Drought is 
expected to impose serious water stress for the growth of MSPs. For the value of REW1-2 m, the 
percentage of severe drought (i.e., REW1-2 m<0.08; Table 2) accounted for 19% in 2013 and 2% in 
2016, but as more as 86% in the first dry year (2014) and 50% in the second dry year (2015). 

Rs was not significantly different among different years, with the mean values varying from 
170.3 to 182.3 W/m? (P=0.40; Fig. 3c). gw in 2013 was 4.88 m but significantly decreased to 5.82 
m in 2016 (Fig. 3d). A slight increase in gw due to heavy precipitation was occasionally observed 
in each year, but overall, gw declined over the entire study period. 

Daily 7;/Ts-max was significantly different (P<0.01) among different years. Daily Ts/Ts-max was 
0.21 in 2013, followed by a decrease of 25% in 2014 and a further decrease of 58% in 2015 and a 
decrease of 52% in 2016. The maximum daily Ts/Ts-max in 2013 was 1.00, which decreased to 0.61, 
0.32 and 0.55 in 2014, 2015 and 2016, respectively (Fig. 3e). 


3.2 REWo-1m change trends in rainless episodes 


A total of 22 rainless episodes with durations ranging from 7 to 37 d were recorded over the 
studied period, in which 8 typical rainless episodes with the distinguishable initial values of 
REW-1 m were plotted to display the change trend of REWo-1 m with cumulative daily ETo (Fig. 
4a). We found that soil moisture decreased with cumulative daily ETo more rapidly for wet soil 
than for dry soil. For instance, when soil was well wetted by heavy precipitation with an initial 
REWo-1 m as high as 0.94 (DOY 197-208 in 2013), the value of REWo-1 m decreased at a rate of 
0.69 per 100 mm of cumulative daily ETo during 12 rainless days. For episodes with an initial 
REWo-1 m of 0.41 (DOY 182-201 in 2016) and only approximately 0.13 (DOY 212-230 in 2015), 
the mean decline rates of REWo-1 m was about 0.25 per 100 mm of cumulative daily ETo during 20 
rainless days and 0.08 per 100 mm of cumulative daily ETo during 19 rainless days, respectively. 
The rate of decrease in the REWo-1 m value with cumulative daily ETo in all episodes were 
negatively related to the decrease in the initial value of REWo-1 m (Fig. 4b). 


3.3 Response of daily T;/ET» to soil water availability 


Daily ratio of T, to ETo (7T;/ETo, dimensionless) was used to elucidate the influence of 
meteorological factors. The T/ETo values during daytime on sunny days varied with soil moisture 
and drought levels were plotted in Figure 5. Generally, T;/ETo annually decreased, and the 
regression slopes relative to REWo-1 m were greater in dry years than in wet years (Fig. 5a). 
Furthermore, mean 7;/ETo annually decreased under the same drought level (Fig. 5b) (P<0.01). 
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Fig. 3 Seasonal dynamics of (a) precipitation, (b) relative extractable soil water (REW), (c) solar radiation (Rs), 
(d) groundwater level (gw) and (e) normalized daily tree water use (daily Ts/Ts-max) during 2013-2016 
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Fig. 4 (a) Relationship between relative extractable soil water in the upper 0-1 m layer (REWo-1 m) and 
cumulative daily reference evapotranspiration (ETo) during 8 typical rainless episodes selected from a total of 22 
rainless episodes during 2013-2016, represented by the day of year (DOY) from (1) DOY 197-208 in 2013, (2) 
DOY 183-193 in 2013, (3) DOY 137-274 in 2014, (4) DOY 246-282 in 2013, (5) DOY 162-198 in 2015, (6) 
DOY 182-201 in 2016, (7) DOY 198-226 in 2014 and (8) DOY 212-230 in 2015; (b) Relationship between 
regression slope of REWo-1 m relative to ETo (from 22 rainless episodes based on linear fitting) and the initial 
value of REWo-1 m. 


The value of T,/ETo ratio under non-drought condition was 0.28 (+0.12) in 2013, and it decreased 
to 0.21 (+0.09) in 2014 and 0.11 (+0.05) in 2015, but slightly recovered to 0.14 (+0.08) in 2016. 
Particularly, after two consecutive dry years, T;/ETo ratio under moderate drought conditions was 
as low as 0.03 (+0.01) in 2016, which was even lower than that under severe in 2015, indicating 
dramatically reduced tree water use after prolonged droughts. Within a year, the daily 7T;/ETo ratio 
at different drought levels showed a significant difference (P<0.01), except in the first dry year 
(P=0.57). 
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Fig. 5 (a) Relationship between daily ratio of tree water use to reference evapotranspiration (7/ETo) and 
relative extractable soil water in the upper 0-1 m soil depth (REWo-1 m) during 2013-2016; (b) T/EToratios at 
different drought levels in each year. Bars indicate standard deviations. Do, non-drought; Dsi, slight drought; Dmoa, 
moderate drought; Dsev, severe drought. 


3.4 Recovery abilities of tree water use after prolonged drought periods 


Two typical sunny days with soil moisture under drought and rewetted conditions after heavy 
rains in each year were used to compare the patterns of hourly 7)/Ts-max. On a slight drought day 
(DOY 282) in 2013, the maximum hourly Ts/Ts-max was 0.55 (Fig. 6a) and slightly recovered to 
0.57 on a non-drought day (DOY 285) after a 42.4 mm precipitation event. In 2014, the maximum 
hourly 7;/Ts-max increased from 0.30 in a moderate drought day (DOY 224) to 0.54 in a 
non-drought day (DOY 238) after an 84.4 mm rainfall event (Fig. 6b). However, the maximum 
hourly Ts/Ts-max On a severe drought day decreased to 0.14 (DOY 228 in 2015; Fig. 6c) and 0.05 
(DOY 196 in 2016; Fig. 6d). Moreover, even if the soil was well rewetted by precipitation in the 
dry years, the maximum hourly Ts/Ts-max in 2015 (DOY 237) and in 2016 (DOY 207) accounted 
for only 44% and 70% of that in 2013, respectively, which indicated a reduced recovery ability of 
MSPs regarding tree water use under prolonged and severe droughts. 
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Fig. 6 Comparison of normalized hourly water use (hourly T:/Ts-max) on a typical sunny (drought) day with that 
on a typical rewetted and drought-free sunny day. It should be noted that these two days selected in a year are 
most likely to be close in time. Squares indicate a drought day (slight drought (Dsi) in 2013, moderate droughts 
(Dmoa) in 2014 and 2016, and severe drought (Dsev) in 2015), and circles indicate a non-drought day after 
precipitation (Do). The arrows indicate the transition from drought to non-drought. 


3.5 Relationship between 7,/ET» and groundwater level 


To further explore the reasons for suppressed tree water use after prolonged and severe droughts, 
we analyzed the relationship between T/ETo and groundwater levels over the studied period. We 
found that 7;/ETo ratio on sunny days generally showed a negative linear relationship with 
groundwater levels (Fig. 7). In the first two years, gw was maintained at a relatively high level, 
varying from 4.88 to 5.27 m, and 7,/ETo ratio varied from 0.17 to 0.38. However, in 2015 and 
2016, gw decreased from 5.32 to 5.78 m, and 7,/ETo ratio decreased from 0.18 to 0.06. The 
suppressed T;/ETo even under well rewetted soil conditions in 2016 was consistent with the 
findings of the continuous decline in the groundwater levels. 
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Fig. 7 Relationship between ratio of tree water use to reference evapotranspiration (T/ETo) on sunny days and 
groundwater (gw) levels during 2013-2016. Bars indicate standard deviations. 


4 Discussion 


4.1 Reduction and incomplete recovery of tree water use after droughts 


In this study, moderate drought in the first dry year (2014) and further severe drought in the 
second dry year (2015) resulted in a marked decline in the water use of MSP. Compared with the 
normal year (2013), the daily Ts/Ts-max was decreased by 32% in 2014 and 62% in 2015. The 
significant reduction in tree water use under drought conditions reflected the conservation of MSP 
in water consumption. This is consistent with the findings on Scots pine trees that have a very 
strong plastic response to intermittent drought (Levitt, 1980; Arneth et al., 2006; Verbeeck et al., 
2007; Mereu et al., 2009). It was reported that the main water source for MSP was from soil 
moisture in the upper soils (Su et al., 2006; Wei et al., 2013; Song et al., 2014). This result was 
supported by our findings that water use of MSP was mainly dependent on soil moisture in the 
upper 0-1 m soil layer. In general, the value of T/ETo ratio linearly decreased with the decline of 
REWo-1 m in each year, and the slopes were more negative in the years within and after severe 
drought than in the normal year (2013) (Fig. 5). Therefore, water availability in the 0-1 m soil 
layer is important for the water use of MSP due to the shallow-rooted traits of MSP that 85% 
roots of MSP are occurred in the 0.0-0.4 m soil layer (Jiang et al., 2002; Zhu et al., 2005; Zhu et 
al., 2008). 

Although rapid recovery of water use of MSP sporadically occurred with increasing soil 
moisture after effective precipitation, the degree of recovery in the dry years (2014 and 2015) was 
lower than that in the normal year (2013) (Fig. 6). However, the maximum daily Ts/Ts-maxin the 
wet year (2016) that was expected to completely recover due to the high annual precipitation 
amount accounted for only 52% of that in the normal year. Water use of MSP in two successive 
dry years could not recovery to the normal year (2013) (Fig. 3e). The incomplete recovery of 
water use of MSP under the non-drought conditions in 2016 could be attributed to the variation in 
hourly 7;/ETo that accounted for 89% of that in 2013. This result is consistent with the reports 
that changes in soil moisture affect water use, but the constraints likely remain over a long period 
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after droughts (Yan et al., 2016). It is previously reported that fine roots of Scots pine die quickly 
under drought conditions (Vanguelova and Kennedy, 2007), which might be the main direct cause 
of the decline in water use capacity even after soil water status was recovered. 


4.2 Effects of groundwater on tree water use under prolonged drought conditions 


A deep taproot enables trees to use water from deep soils and/or groundwater under drought 
conditions (Barbeta et al., 2015; Hentschel et al., 2016). MSP is a dimorphic-rooted species with a 
maximum taproot depth of 5.20 m, which was observed in a 42-year-old tree in sandy soil near 
our site (Jiang et al., 2002). The water absorption by the deeper taproots of MSP was measured by 
isotopic technique when soil water content in the upper 0.6 m soil layer declined to 3.6% (Wei et 
al., 2013; Song et al., 2016a). This finding was supported by our results of the extreme depletion 
of soil water in the 1—2 m soil layer. Meanwhile, gw value decreased from 4.88 to 5.82 m over the 
studied period, and in 2014, gw value was much lower than the maximum depth of taproot (5.20 m). 
Considering the findings that 7;/ETo ratio during the daytime on sunny days decreased linearly 
with the decline in groundwater levels (Fig. 7), we concluded that significant decline in 
groundwater levels was partly responsible for the reduced tree water use even when soil water 
availability was gradually recovered. The impact of groundwater on tree water use was most 
prominent under extensive and prolonged drought conditions, which was greatly but incompletely 
replaced by rewetted upper soil. The results from MSP in this study are consistent with the report 
that dependence of trees on groundwater increases as precipitation decreases (Kume et al., 2007). 


4.3 Effects of soil texture on soil drought 


Sandy soils are characterized by a low water-holding capacity and high hydraulic conductivity. 
Therefore, these soils provide a harsh environment for plant growth. Although soil moisture in the 
upper layer was easily recharged from the precipitation, it rapidly decreased in a short term even 
at higher initial soil moisture conditions in the sandy soil due to the soil texture and strong 
evaporation (Fig. 4). Therefore, droughts tended to occur more frequently in the sandy soil due to 
the high inter-monthly variations in precipitation. Moreover, coarse sandy soil effectively 
prevents capillary rise of water (less than 0.5 m) (Vincke and Thiry, 2008). It is reasonable that 
the contribution of groundwater to upper soil moisture by capillary rise is non-significant due 
to >91% sand particles in the upper 2 m soil layer in the studied plot. As a result, more than half 
of days in the dry years (2014 and 2015) were drought, which was believed to induce persistent 
water stress on the growth of MSP. 


4.4 Implications for management of MSP plantations under climate change 


T;/ETo ratio reflects the response of tree water use to atmospheric water demand, which allows us 
to focus on determining the effects of soil water availability on tree transpiration. In our study, 
T,/ETo ratio generally declined with a decrease in the upper soil moisture. T;/ETo ratio also 
showed a significant difference between different drought levels and recovered to a certain extent 
when the whole soil layer was thoroughly rewetted in 2016 (Fig. 5b). This study reported the 
plasticity of MSP under severe drought conditions through the drastic reduction in water use. 
Additionally, the lateral roots of an MSP tree can extend to lengths that are approximately 0.65 
times that of tree height (Jiang et al., 2002; Su et al., 2006), which also helps MSP trees to 
efficiently obtain water from the upper soil layer (Song et al., 2014). However, incomplete 
recovery of water use revealed that prolonged and severe droughts might greatly limit the ability 
of trees to adapt to climate change. In addition to random drought episodes that can greatly reduce 
water use, the progressive decline rate in the groundwater level (1 m/a) induced by extensive 
and/or intensive agricultural land use and broadleaved forests plantation in the study area (Song et 
al., 2016b) also prevented the recovery of water use, both of which would cause a deficit in the 
water supply and serious degradation of forest stands if allowed to continue for two or more 
growing seasons. As a countermeasure against the ongoing degradation of MSP, it is more 
effective and more urgent to prevent unrestrained exploitation of groundwater by reducing 
agricultural land or converting into other land use such as coniferous forest and shrubland than by 
merely adjusting stand density of MSP. 


5 Conclusions 


High inter- and intra-annual variations in precipitation and soil moisture in sandy soils led to 
frequent droughts in the semi-arid Horqin Sandy Land. Consequently, water use of MSP strongly 
responded to soil water availability, especially in the upper 1 m soil layer, as these trees are 
shallow-rooted species. Severe and prolonged soil drought greatly decreased water use of MSP in 
a short term. The progressive decline in groundwater restricted the degree of water use recovery 
in MSP even if the soil was rewetted by precipitation. Our results suggested that the degradation 
of MSP plantations is attributable to the high temporal variations in precipitation, low soil water 
retention and declining groundwater level. In addition to adjusting planting density and stand 
structure, preventing groundwater level from decreasing by reducing irrigated farmland area is 
recommended as an alternative method to stop ongoing degradation of MSP in the Horqin Sandy 
Land of northern China. 
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